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Lagrangian  Multipliers  and  Superfluous  Variables 


1.  Introduction 

This  paper  uses  superfluous  coordinates  and  the  "constraint  jiggling" 
approach  (usually  introduced  when  performing  displacements  violating  con¬ 
straints  in  order  to  compute  the  forces  of  constraint)  to  develop  a  new  form  for 
Lagrangian  Multipliers.  As  a  part  of  the  development,  this  paper  generates 
unified  equations  combining  a  new  variant  of  the  Euler-Lagrange  equations  with 
a  new  closed  form  expression  for  these  multipliers  in  terms  of  the  superfluous 
coordinates  and  the  boundary  conditions.  These  results  elucidate  the  mechanism 
by  which  the  constraints  determine  the  Lagrangian  Multiplier,  and  their  role  in 
the  Euler-Lagrange  equations  with  constraints.  This  approach  facilitates  the 
computations  for  the  following  reasons: 

(a)  Finding  the  inverse  functions  requires  the  solution  (at  worst)  of  a  system 

of  m  equations  in  3n  coordinates,  where  the  3n-m  independent  coordinates 
are  treated  as  known  constants  and  no  derivatives  are  involved.  Once 
done,  the  dependent  coordinates  and  their  derivatives  can  be 

replaced  in  the  constrained  form  of  the  Euler-Lagrange  equations  derived 
in  the  paper  by  ...  ,gm  and  their  derivatives,  so  these  equations  will  now 
contain  only  the  independent  )n~m  coordinates  and  their  derivatives.  The 
standard  treatment  requires  the  solution  of  a  system  of  3»  differential  equa¬ 
tions  and  m  algebraic  equations  in  3«+m  variables. 

(b)  Once  these  Euler-Lagrange  equations  are  solved,  the  Lagrange  Multipliers 
are  obtainable  without  further  solution  of  equations. 

(c)  The  inversion  can  often  be  performed  by  inspection,  as  is  the  case  in  the 
example  cited  in  the  paper. 

(d)  A  good  choice  of  coordinates  may  trivialize  some  of  the  m  equations  when 
computing  the  inverse  functions. 

2.  Statement  of  the  Problem 

We  assume  that  there  are  *  particles  traveling  in  3-dimensional  space,  sub¬ 
ject  to  m  smooth  independent  holonomic  constraints 

•••. **•;<)-«.  for  a  -  1 . m  (l) 

We  treat  this  problem  in  the  3n-dimensional  configuration  space  of  the  n  parti¬ 
cle  system,  by  concatenating  the  *  3-dimensional  position  vectors  into  one  vec¬ 
tor  in  Ru.  Since  the  constraints  are  assumed  to  be  smooth,  the  constraint  force 
corresponding  to  constraint  surface 
f ei*\i  ■  ■  ■  p***;*)™*® 
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most  be  co- linear  with 

«c*i.  •  •  •  •  ■  •  >i3«iO  Ui 

where  4  is  the  k*  canonical  unit  vector  in  Ru.  Therefore,  the  total  constraint  force  & 
most  be  a  linear  combination  of  the  tya  . 

We  also  assume  that  the  constraint  forces  (  and  gradients  )  are  locally  linearly 
independent  along  the  trajectory  of  the  system  in  configuration  space. 

Note  that  all  vectors  in  this  paper  that  do  not  involve  time  will  be  Ru  dimen¬ 
sional,  whereas  vectors  that  do  involve  time  will  be  A3"*1  dimensional. 

3.  Results  to  be  proved 

The  following  discussion  will: 

(1)  Derive  a  new  expression  for  the  Lagrangian  multipliers 

\  r*  i  f  dL  ^9a . 

#u  a  *[0ioj~3*a  »t 

where  the  fJth,  . . .  . . .  ,**,;»)  are  the  inverse  functions 

of  the  m  constraint  functions  /^z,, . . .  where 

*  /«(* i,  •  •  •  .  for  a  -  l,  .  .  ,  m. 

(2)  Derive  a  new  variant  of  Lagrange’s  equations,  using 
the  superfluous  coordinates. 

(3)  Show  that  the  \f(t)  are  components  of  a  generalized 
constraint  force  in  R3m. 


4.  Finding  the  Dependent  Coordinates  and  Relating  them  to  the 
Independent  Coordinates 

Let  (s„  .  .  .  ,xuX)  be  an  arbitrary  point  of  Ru+l.  We  next  want  to  find  the 
dependent  coordinates,  and  show  that  we  can  relate  them  (  in  a  neighborhood 
of  (x„  .  .  .  ,*„;!)  )  to  the  independent  coordinates.  The  dependent  coordinates 
are  also  referred  to  as  "  superfluous  "  in  some  of  the  physics  literature.  We 
have  assumed  that  the  m  constraints  are  independent,  and  therefore  have 
assured  that  there  are  3n  -  m  independent  coordinates.  We  will  intentionally 
remain  vague  about  the  exact  meaning  of  "neighborhood"  in  R3n,  but  a  good  dis¬ 
cussion  may  be  found  in  the  first  cited  reference. 

Note  that  while  it  is  usually  possible  (  discussion  of  "usually"  deferred  )  to 
relate  the  dependent  to  the  independent  variables,  this  can  only  be  done  locally 
(  therefore  the  "neighborhood"  ).  Variables  that  are  dependent  in  one  area 
could  in  principle  be  independent  elsewhere.  In  any  case,  the  functional  expres¬ 
sions  relating  the  dependent  to  the  independent  variables  may  only  have  local 
validity.  In  the  examples,  it  is  shown  that  in  some  cases,  this  locality  is  not 
much  of  a  problem,  because  the  neighborhoods  are  really  enormous.  The  paper 
will  generally  assume  some  neighborhood  N  in  which  the  discussion  takes  place. 
The  Inverse  Function  Theorem1  provides  the  necessary  conditions  for  the 
existence  of  the  m  dependency  relations  (  the  "usually"  mentioned  above  )  and 
for  the  existence  of  neighborhood  N. 


We  start  by  defining  a  function  that  is  a  set  of  3n  coordinate  transformations  with 
tune  as  an  added  coordinate.  We  next  show  that  a  certain  determinant  ( its  Jacobian  ) 

is  non-zero  in  a  neighborhood  of  (*, . This  condition  is  required  by  the 

Inverse  Function  Theorem.  This  theorem  will  then  justify  the  existence  of  the  inverse 
coordinate  transformations  and  the  appropriate  neighborhood  N. 

We  pick  a  point  Cct . and  we  assume  that  the  m  constraint  functions 

/■Ui. . . .  .jtj „u)  for  a  -  1. ...  .m  are  independent,  i.e.  have  gradients  Vfx . Vfm 

which  are  linearly  independent2  at  the  point.  Hence  the  matrix 


ill.  ILl 

d*x  dtu 

9/m  9/m 

9*  i  9t 

has  m  independent  columns  (  or  equivalently3  has  rank  m  ). 

By  relabeling  the  variables  if  necessary,  we  may  assume  that  the  leftmost  m 
columns  are  linearly  independent,  so  that 


where 


/,(*»,  •  •  •  ,*s.;0  ■  */.  for  j  -  m+l . 3». 


The  Jacobian  of  F  is  the  (3*+l)  by  (s*+i)  determinant 


dx, 

9fi 
9x |a 

dt 

ill 

ill 

9fu 

9fu 

—  det 

9*i 

9*» 

9*i 

9tu 

dt 

9fm 

9fm 

0 

0 

1 

9*i 

9*m 

Therefore  the  Jacobian  of  F  is  non-zero  in  the  aforementioned  neighborhood  of 

(Sl>  •  •  ■ 

Define  jr«  -  fj.*\,  •  ■  •  ,  *s,;0>  for  a  —  1,  ....  in. 

The  Inverse  Function  Theorem  now  allows  us  to  conclude  that  there  are  3n+i 
functions 

•  •  •  #  ■  ■  ■  ,  *U,0  for  •  -  1,  ...  ,  3l»+l 

1)  Possessing  continuous  partial  derivatives. 

2)  Defined  on  points  (f„  .  . .  in  F[N),  the  image  of  some 

neighborhood  N  of  (*,,  .  . .  ,xU'D- 

In  addition,  these  functions  satisfy  the  important  (  inverse  )  condition 

#»(/ l{*U  •  •  •  •  >  f m(*l>  •  •  •  •  •  ■  >*3ti0  “  **, 

for  t  ■  l,  .  .  .  ,3«+l,  where  x,.*.,  -  t. 


These  functions  are  the  required  -inverse  transformations,  and  will  be  used 
extensively  in  nhe  following  sections. 

Recall  that  /«(*,,  .  .  .  , xs.;f)  •  for  a  -  i . m,  due  to  the  constraints. 

Therefore 

*»  ■  fc(eli  •  •  •  i  em*m+li  •  •  •  » 


SO 


for  K«,  xt  depends  on  xw+,,  .  .  .  ,zs. 

for  k>m,  ft  is  the  projection  whose  value  is  x*' 


These  observations  are  crucial  to  the  rest  of  the  argument. 

Note  that  since  the  t»  functions  associated  with  the  independent  coordi¬ 
nates  are  projections,  only  the  m  functions  associated  with  the  dependent  coor¬ 
dinates  must  be  found.  Note  also  that  these  functions  may  sometimes  be  found 
by  inspection,  as  may  be  seen  in  the  example  cited  later  in  the  text.  Even  when 
these  functions  must  be  derived  by  solving  equations,  the  3 n-m  independent 
coordinates  may  be  treated  as  known  constants,  so  only  an  m  x  m  system  must 
be  solved.  Lastly,  a  good  choice  of  coordinates  may  trivialize  some  of  the 
remaining  m  equations. 


5.  Path  Variation  and  Coordinate  Independence 

The  next  argument  depends  on  Hamilton’s  Principle4.  We  will  use 
Hamilton's  Principle  in  N  to  obtain  generalized  Euler-Lagrange  equations  for 
the  *  particle  system  with  the  stated  constraints. 


4 


Recall  that 


*«  -  *«(*»+!#  •  ••»*«■;*)  for  or  -  1,  ....  m , 
so  are  functions  of  the  independent  coordinates  *m+1) 

We  define 

•  •  •  f  **•  •  •  ■  >  *S»»0  ■  •  •  •  i  9mizm+V  •  •  •  >XS*’9\>  ■  •  •  •  ■  ■  >*3*>0' 


Also  note  that  &i„  -  2  TZmSxn  and  that  the  path  variations  in  Hamilton’s 

y— w  9xi 

Principle  satisfy  the  condition 


&*j  -  |d*yj  for  j  -  m+1 . 3« 

Analogously,  we  define  5 ja  m  so  we  can  compute  that 


We  now  use  Hamilton’s  Principle  on  L„  for  points  (*.+lI  .  .  .  ,**»;*)»  where 
(*i,  •  ■  • ,  *u',t)  is  in  neighborhood  N. 


0  m  Sfdt  f«i(  *«,+!,  .  .  .  >  *S«|*«I+Ii  •  •  •  >  *3«i0 


r  ..  d  L  -  d L  OL  c  o  JL  r  • 

-Jdt  Zj-z S7rdz>  +  S7 St—5*, 

«l  y-m-M  <'*y  I»i5la  y-w-H  ^*y 


r,.  1  A  dL  M,  ^  dL  ,  ^  dL  </  I, 

/“  1  Zijr  ^  ir4*/  +  ^  aT6z/  +  Stt-  £  tt  t— 


*  tf-  £  | 'f‘i<  *  fi  f-H, 

*-l  u*i  ;-m+l  0*; 


-J*  E  E- 


dL 

d 

99c 

dx„ 

il 

,  dxi  . 

.  M  Aaz, 

+  £  £  a"  +  —  6xi 

y-*-w  l»-l  d*„  #*/  d*y 


I'll  I  If-WW, 'I  H..L 


We  now  recall  the  following  observations: 
for  j  -  m+1, . . .  ,3» 

2)  «.,(*,)  -0 

3)  $*/0a)  “  ®* 

We  use  the  above  observations  while  integrating  the  dx,-  term  by  parts  in  the 
equation  preceding  the  observations  to  obtain: 


*« 

0-  Jit 


3m 

E 

jmm+l 


A  dL 
“J  dx.  dx,- 


.  Ik 

I**,-. 

9*j 

&Xj 


u 

-  E 

/—+« 


hit 


dL 

9ga 

+  -f 

2k 

Bxi 

dx  a 

d*i  hi  dx.  it 

,  dxi  . 

dt 

dx  j 

3 

Note  that  the  term  £•—  - 


-i  di„  it 


dga 


dx, 


•)  ) 


**  Sa 

o-fit  £ 

it  y— i+i 


cancels  out,  and  we  now  have  the  equality 

k 


m 

dL  i 

dL 

dL  i 

dL  ' 

L 

dx.“  it 

dx, 

dx,-  it 

.  dii 

Now  the  Si,  for  /  —  m+l.  .  ,  . ,3«  are  independent  (  while  the  6x .  ix, 


dL 

dXj 


(  dL 

m 

+  V1 

r  » 

dL  i  [  dL 

\9ii, 

+  Li 
«■! 

it  [  dx. 

.  . 

dfa 

dXj 


«  0  for  j  m  m  +  1,  .  .  .  ,3n. 


We  will  now  use  this  equation  to  motivate  our  formula  for  the  Lagrangian  Mul¬ 
tipliers  and  to  derive  the  usual  expression  for  the  Lagrangian  Equations  with 
constraints. 


6.  Deriving  the  Expression  for  the  Lagrangian  Multiplier 
Let  us  rewrite  the  previous  equations  in  the  form 


i 

’  dL  1  dL  A 

— 

dL  | 

dL 

it 

**  r  **  i 

dt  j 

J 

- 

for  j  -  m+1,  . .  .  ,3n  which  is  valid  for  points  (*„  .  .  .  in  N. 
Now  let 

•  •  •  >  *hi0 


®  #«^/ l(*l»  •  •  •  i  *s»>0«  •  •  ■ 

•  xm  for  a  «■  1,  .  .  .  ,  m 

./»(* i,  • 

•  •  »  •  * 

•  >  *S»i t) 

Therefore  -  0,  for  j  - 

dx, 

m+1,  .  . 

•  t  3»  9 

since  no  change  in  x,-  affects  the  value  of  G*. 

Also,  recall  that  jr.  -  /*(* ,,  .  .  for  a  «  i,  .  .  ,m. 


(2) 


6 


It  follows  that 


or 


dG* 

dXj 

dg* 


A  dg* 

,  d3<*  dh 


dg* 

dx> 


9*i  9*i 

Combining  this  result  with  equation  (2),  we  see  that 


± 

it 


dL 


d*i 


dL 

9xj 


-E 

#■1 


i_ 

it 


dL 


dia 


dL 

dx„ 


for  j  —  m+ 1, 


,  3n 


dja 

dxi 


m 

„  v 

i 

dL  | 

l 

dL 

j 

A  dga  df0 

~  1 

it 

dxa) 

~dx* 

.  h  d*9  dxi  . 

'  E 


Ill 

± 

dL 

1  dL 

dg* 

E 

it 

dxa 

1  d*a 

dgp 

df9 

dx. 


(3) 


So  if  we  define 


MO  ■  E 


dL 

dL 

it 

diaj 

dr* 

dg* 


]d1t 

we  note  that  \f(t)  is  the  component  of  the  generalized  force 

j  nr 

9i,  ■  ■  ■  ,vm  corresponding  to  the  m-dimensional  force  —  - 

it  dx* 

z,,  ...  ,xm  coordinates.  Recail  that  x,,  . 


(4) 


in  coordinates 

in  the 
dr* 

,xm  are  the  dependent  ••oorciinates  i.e. 

the  coordinates  that  are  functions  of  xm+1,  .  .  .  ,xu  due  to  the  constraints). 

We  now  summarize  the  results  of  the  last  two  sections  by  combining  the 
equations  for  the  Euler-Lagrange  equations  and  the  multiplier  equations: 

(MO  for  1  m  l.  • 


3« 

E 


i 

'  dL 

dL 

dg »  , 

it 

dr„ 

dyi  ” 

for  I  —  m+1. 


,  m 
.  ,  3n 


(5) 


7.  Interpreting  the  Lagrangian  Multiplier 


We  can  interpret  the  meaning  of  — 

at 


dL 


dxk 


dL 

dx„ 


by  appealing  to  a  basic  result  in 


Lagrangian  Mechanics,  if  we  assume  smooth  constraints  and  conservative  forces. 

Recall  that  if  C*f-i>»2.  Qxi-w  are  the  components  of  the  total  force  on 

particle  q,  then  £  is  the  total  force  in  configuration  space,  and  we  have 

d_  87  dT  m  Q 

it  8ij  dxk  *' 


for  k  *  1 . 3/i, 


(6) 


where  T  is  the  kinetic  energy. 


Now  <?*  •  F*  +  /?*,  where  Fk  is  the  ka  component  of  the  3n-dimensional  external 
force  and  Rk  the  ku  component  of  the  3n-dimensional  smooth  constraint  force. 
We  assumed  conservative  forces,  therefore  there  is  a  potential  V  such  that 

Ft  9*k  * 

We  proceed  to  define  LmT-V  as  usual,  so  by  (6), 

-f  (— “  R»  3». 

it  9xk 

Utilizing  this  result  for  the  first  m  variables  only,  we  see  that 

MO  - 

<wl 

so  the  Lagrangian  Multiplier  appears  to  be  the  portion  of  the  total  generalized 
constraint  force  due  to  the  dependent  coordinates.  However,  since  the  con¬ 
straints  are  smooth,  the  constraint  force  has  zero  components  in  the  directions 
corresponding  to  the  independent  coordinates,  and  therefore  the  components 
corresponding  to  the  dependent  coordinates  in  fact  determine  the  total  force  of 
constraint.  Note  that  the  equations  (5)  combining  the  multipliers  and  the 
Euler-Lagrange  equations  have  the  constraint  force  on  the  right  hand  side. 

8.  Showing  Multipliers  Satisfy  the  Usual  Constrained  Euler-Lagrange 
Equations 

Now  for  j  -  m+l,  .  .  ,3n  and  points  in  N,  we  immediately  have  by  (3)  and  the 
definition  of 
i  9L  9L  f,.  ,.Jfs 

-°- 

Now  for  7  »  l,  ....  m: 

Again,  for  a  -  l, . . . ,  m 

9cHf  i(*i»  •  •  •  •  •  •  >/• i(*i>  •  ■  ■  •  *s» ■  •  >  *s*»0  “ 

so  if  -y  #  o,  *7  does  not  affect  the  value 
of  to.  It  follows  that: 

f 'ilz-Le.  _ 

d*i  “  d*i 

-  the  Kronecker  delta. 

Therefore  we  see  that 


4  TT  -TT-  tW- 


it  [  9i 7  J 

_  J-\lL.)±L  _  i  \  9L)  dL 

it  [  dxn  J  &*■,  it  9x1  1 

Combining  the  results  for  j  and  t,  we  have 

EMO--  -0,  for  i  —  1 . 3n . 

This  is  the  usual  form  of  the  constrained  Euler-Lagrange  equations, 
tt.  Relating  the  Forces  of  Constraint  and  the  Lagrangian  Multipliers 
We  showed  that  —  (•—*] — —  -  Rk,  the  ***  component  of  the  constraint 


it  [  9ik  J  9xk 

force  when  the  external  forces  are  conservative  and  when  the  constraints  are 


smooth.  Therefore,  Rk 


d_ 

dt 


'  dL  |  dL 
dxk  dxk 


£wo 

^-i 


dxk  ' 


We  can  rewrite  this  in  3 n 


dimensional  space  as  it  =  £Xf»(0^/p>  which  coincides  with  another  popular 
expression  for  the  constraint  force. 


This  expression  makes  physical  sense  because  smooth  constraints 
must  be  co-linear  with  the  gradients  of  the  constraint  equations. 


10.  Example  for  Constrained  Time-Independent  Two-Dimensional 
Motion 

Now  for  an  example  of  single  particle  motion  in  two  dimensions.  Let  a 
bead  slide  frictionlessly  along  a  wire  whose  shape  is  given  by  the  equation  y=x2. 

Let  the  constraint  equation  ( in  the  article’s  notation  )  be 

o  *  «i  -  /(* i,*2)  -  *2-*?- 


Define 


F(xlfX2) 


*2-*? 

,  *2  . 


The  Jacobian  of  F  is  2*„  so  the  Inverse  Function  Theorem  applies  when  x,*o. 
The  Inverse  Function  Theorem  only  guarantees  the  existence  of  an  inverse  func¬ 
tion  but  doesn’t  provide  its  form.  We  need  a  function  g  satisfying 

g(f{xx,x 2),x2)  -  *t.  This  is  equivalent  to  finding  a  function  g  satifying 
j(x2— x?,x2)  »  xt.  Note  that 

•jgn(x2— (x2— xr))"v~xo— ixn— x7)]  =  x.  if  x,  >  ■) 

»yn((x2— xf )— x2)"\/ 1  Xo— (xo— xf)|  =  xl  if  x,  <  0 


where 


sgn(x) 


1  if  x>0 
0  if  x=0 
•1  if  x<0 


So  we  define  two  functions 

9 1(?.*2)  *  Jgn(x2-y)Vjx2-y  ! 

?2<».*2)  *  »gn(y-x 2)V|x2-y| 

where  y  m  /(x„x2). 

The  reader  can  easily  show  that  if  x,  >  0  then  Ji(/(*i,x2),x2)  =  x,. 

We  define  the  neighborhood  N,  by  Nj  *  {(x,,x2)|  i!  >  0 },  where  is  defined  on 
F[N,]  *  {(/(*i,*2),*2)!(*i.*2)«Ni}-  Note  also  that  is  differentiable  on  F[N,|. 
Correspondingly,  if  x,  <  0  then  g^fix^)^)  -  xx. 

We  define  the  neighborhood  N2  by  N2  s  {(x,,x2) J  xt  <  0  },  where  g2  is  defined  on 
F[N2j  *  {(/ (x,,x2),x2)  1  (xi,x2) t N2}.  Again,  j2  is  appropriately  differentiable  on  F[N2j. 

Note  that  both  of  these  neighborhoods  are  in  fact  quite  large.  In  the  rest 
of  this  section,  we  will  use  g  interchangeably  for  either  gi  or  g2,  depending  on 
whether  (y,*2)<F(N,]  or  (y,x2)«F(N2]. 
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We  see  that 


L 


ym(i,2+if]-myx2. 


The  reader  can  show  that  since  x{  is  the  dependent  variable 

_d_ 
it 

and 


dL. 


dL 

dx i 


mi. 


\ 


it 


M 

diy 


dL  dg  , 
«*! 


mxl 

2*i 


This  is  the  same  value  for  X  as  would  be  obtained  from  the  usual  constrained 
Euler-Lagrange  equation 


it 


dL  ]_  dL  _xdf 
dil  dij  3i, 


0 


with  the  constraint  equation  x2-t?  -  0. 

We  see  that  for  xt  =»  0,  X  is  undefined,  as  should  have  been  expected  from 
the  necessity  for  the  exclusion  of  such  points  in  the  definition  of  N,  and  N2. 


11.  Example  of  Constrained  Time-Dependent  Two  Dimensional 
Motion  (No  Gravity) 

Let  an  infinite  rod  with  a  pivoted  end  at  the  origin  of  a  plane  be  rotating 
in  the  plane  with  frequency  w  about  the  origin.  Let  a  bead  slide  frictionlesslv 
on  that  rod. 

The  constraint  equations  are 
i— rcoe  ut  w  0  m  ex 
y— rain  wt  ■»  0  ■  e2. 

We  translate  to  the  article's  notation  as  follows: 

V\  »  /i(*i,*s»*s;0  -  *i-*3cos  u t 
*  !A* i.*2.*s;0  -  *2-*jsin  w t. 

Note  that  i,  and  i2  may  be  chosen  the  dependent  variables,  upon  examination 
of  the  Jacobian. 


We  next  define 

f/i(*i.*j.*3;0 


/  s(*i.*j»*3;0 

*3 


u 

Xx—  I3COS  w t 


xy—  ijsin  wt 
*3 


U  J 

The  Jacobian  of  F  is  equal  to  1  ,  for  all  points  (i,,i2,i3;f),  so 
we  may  choose  N  m  {(i,,i2,i3;t)}. 


We  need  to  find  functions  9\{yx,y2,x3;t)  and  y^yi.i^aiO  satisfying 
y»(/i(*i.*9.*3;*)./s(*i.*a«*s;0.*s;0  “  »i(*i-*3c°®  wt^-x^m  wf,*3;0  *  *t  and 
fs(/i(*i.*8>**;0./2(*it*2»*»;0»*8;0  *  ?2(*i-*3c°s  wf,z2-x3 sin  ut,x3lt)  =  x2. 

The  obvious  choices  are 
#i(yi»>i.*s;0  -  jri+*sc<*  and 
ifauWai*)  m  ya+xjsin  wf. 

It  is  easily  seen  that  yv  and  y2  are  appropriately  differentiable 
and  that  the  required  conditions 
9i(.fi(*i,**,*»;t),fd*i>*s>x3>t),x3;t)  -  *1 
Jij(/l(*l»*2j*3»0>/2(*li*2>*3>0>*3i0  “  *2* 

are  satisfied. 

Since  L  »  -|-m(ii2+x|],  we  see  that 


_J_ 

it 


dL  I  dL  t 
dik  j  9xk 
for  Jr  —  1,2. 
Consequently, 


mxk 


x,(0  - 

_J_ 

'  dL  ' 

dL 

it 

dxx 

!  » 

j 

dL 

A 

dL 

l  ^ 

jt 

di2 

dx2 

dgl 


dgx 

dg2 


ay. 


“  m*i!  »i-« traps' 

Using  the  fact  that  yt  »  e{  m  o,  we  see  that 
x,  =  r3cos  wf.  We  may  use  this  fact  to  compute  s.  and  obtain 
Xi(f)  *  mx3cos  Jl-2mwijsin  wt  — mw’xj cos  -Jt. 

Analogously, 


X2(t) 


i_ 

dL 

dL 

d9i  , 

it 

dix 

dxx 

3y2  <*!«.. »X-<3 

Y 

_dL_ 

0?2  , 

jt 

3i2 

dx2 

dy2  '*»-*  *»*-*«* 

Proceeding  as  before,  we  see  that 


X2(t)  -  mar's! 


'bfa-«a* 


so  since 


y2m  c2m  o,  we  have  that 

z2  ■»  Xjsin  wt. 

We  can  now  show  that 
X2(f)  —  mx3+2mwx3cos  wt 
— mu^xjsin  w f. 

Now  the  usual  constrained  Euler-Lagrange  equations  are: 

d  ' 


it 


jo.  _iL.Xl,l)f£i.M„»n 

ai,  9«,  'a,,  a,. 


and 

it  {  di2 
with  constraints 
0  >■  xt— x3cos  wt  and 
0  »  x2— Xjjsin  w t. 
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A-l. 


The  reader  can  show  that  we  get  the  same  equations  for  x,(<)  and  x2(<),  so  the 
two  methods  are  equivalent. 

Note  that  the  functions  gug2  were  obtained  by  inspection,  so  that  the  mul¬ 
tiplier  values  were  obtained  essentially  without  solving  equations. 

12.  Example  Comparing  the  Standard  and  New  Approaches 

By  way  of  illustration,  let’s  compare  the  solution  processes  required  to 
obtain  the  equations  of  motion  and  forces  of  constraint  for  a  particle  of  mass  m 
moving  on  a  surface  defined  by  the  constraints 

A(*i»*a»*s)  ■  *i*2*s  -  «i  9*  o 

and 


subject  to  the  potential 


/ *  *i+*a  -  «2 
V(*l.*2.*j)  -  J*(*l2+*22+*i2)- 


We  define  new  coordinates  by  gx  -  /i(*1(x2,*3)  and  y2  *  /s(* v,x2,x3).  The 

Lagrangian  now  has  the  form 

i>(*i.*2.*s.*i.*2,*j)  -  jfc(*i2+*22+X,2). 


It  may  easily  be  shown  that 


-i 

dl  ] 

_  dl 

•it 

,)z ,  j 

0zx 

dl 

dl 

dt 

dz2 

jL 

'  dl 

dL 

dt 

. di 

dx. 

mi",+£x. 

mx2+fcr2 

mij+bj 


The  standard  formulation  provides  the  equations 


0 

0 


—  ^mx^+Jfc*!  j— Xjx«xj—  X2 

—  ^mx2-)-itx2j— X,x,x3— X2 
0  -  |mx3+/fcx3j-X,Z|*2 


x,x2xs  -  ex 
Xi+x2  “  e2, 

which  requires  solving  5  equations  in  5  unknowns. 

The  formulation  in  this  paper  requires  the  computation  of  the  inverse  func¬ 
tions  in  (say)  N,  m  {(xlIx2,xJ)jxl>x2},  the  case  for  N2  *  {(x,,x2,x3)ix,<x2}  is  identical 
except  for  interchanging  the  definitions  of  the  functions  described  below. 

The  inverse  functions,  obtained  by  solving  two  equations  in  two  unknowns  are: 

?2+  [»22-4yi/x3j2 
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-  - 2 - • 

Note  that  x,  -  ji(jfi,»2,*s)  ,  *2  -  »2(ifi,ifj.*s)»  and 

*s  -  ?j(»l.»**s)- 

We  now  compute  all  first  order  partial  derivatives  of  gt,gt,  and  gs  evaluated 
at  (tfuej.xj)  and  substitute  these  values  in  the  new  equations,  immediately 
obtaining  the  following  system  of  equations: 

(  ..  1  [ca2— 4ei/*s]  *  /  ",  [e22-4e,/zj]  2 

«■  —  1  mxl+kxl  +  I  mij+ixj  - 

'  '  *j  v 

\i- j(m*,+fcr,)[1+c2[esa-4e1/x,]  2  ]+y  (mij+b,)  [l-e2  [e22-4c,/*s]  2] 


[mij+fat,] 


«1 

«22-4<i/*s 

1 

2 

e22-4e,/xs 

l  *>2 

mxj+ix  j  1 

J 

+«is+b3. 


We  see  that  this  approach  involved  finding  two  inverse  functions  by  solving 
a  2  by  2  system  of  equations,  and  then  computing  somewhat  harder  derivatives 
than  was  required  in  the  standard  case.  On  the  other  hand,  these  equations 
contain  only  r,  and  its  derivatives  once  x,  and  j?  and  their  derivatives  are 
replaced  by  ji(cx,c«,x3)  ana  j^cx,co,x.,)  and  their  ierivatives.  in  addition.  •  nt- 
Lagrange  Multipliers  were  available  immediately.  Therefore,  we  see  that  solving 
the  equations  to  obtain  inverses  and  substitution  of  the  inverse  functions  and 
their  derivatives  for  the  dependent  coordinates  and  their  derivatives  replaces 
solving  a  more  complex  system  of  equations. 


13.  Conclusion 

This  paper  has  related  the  Lagrangian  multipliers  to  superfluous  coordi¬ 
nates,  and  has  derived  a  new  equation  combining  a  variant  of  the  Euler- 
Lagrange  equation  with  a  new  closed  form  for  the  multipliers  in  the  process. 
Since  the  new  multiplier  equation  lead  to  the  usual  form  of  the  Euler-Lagrange 
equations,  the  new  form  for  the  multiplier  must  indeed  yield  the  same  value  as 
the  usual  Lagrangian  multiplier. 

The  author  feels  that  this  approach  improves  one’s  understanding  of  the 
meaning  of  the  multipliers  and  the  role  of  (and  relationship  to)  superfluous  coor¬ 
dinates.  In  addition,  this  approach  facilitates  the  computation  of  the  Lagrange 
Multipliers  by  reducing  the  complexity  of  the  equations  that  must  be  solved,  as 
was  the  case  in  the  example  discussed. 
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